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ABSTRACT 

This article presents a multi-wavelength study towards the massive star- 
forming region S87, based on a dataset of submillimeter/far-/mid-infrared (sub- 
mm/FIR/MIR) images and molecular line maps. The sub-mm continuum emis- 
sion measured with JCMT/SCUBA reveals three individual clumps, namely, 
SMM1, SMM2, and SMM3. The MIR/FIR images obtained by the Spitzer 
Space Telescope indicate that both SMM 1 and SMM 3 harbor point sources. 
The J = 1-0 transitions of CO, 13 CO, C 18 0, and HCO+, measured with the 
13.7m telescope of the Purple Mountain Observatory, exhibit asymmetric line 
profiles. Our analysis of spectral energy distributions (SEDs) shows that all of 
the three sub-mm clumps are massive (110 — 210 M ), with average dust tem- 
peratures in the range ~ 20 — 40 K. A multi-wavelength comparison convinces 
us that the asymmetric profiles of molecular lines should result from two clouds 
at slightly different velocities, and it further confirms that the star-forming ac- 
tivity in SMM 1 is stimulated by a cloud-cloud collision. The stellar contents 
and SEDs suggest that SMM 1 and SMM 3 are high-mass and intermediate-mass 
star-forming sites respectively. However, SMM 2 has no counterpart downwards 
70 /mi, which is likely to be a cold high-mass starless core. These results, as 
mentioned above, expose multiple phases of star formation in S87. 

Subject headings: submillimeter — infrared: ISM — ISM: individual (S87) - 
ISM: molecules — stars: formation 

department of Astronomy, Peking University, Beijing 100871, China 

2 Current address: National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, 
China 



-2- 



INTRODUCTION 



Massive stars play an important role in the evolution of the interstellar medium (ISM) 
and galaxies; nevertheless their formation process is still poorly understood from the ob- 
servational perspective because of their relatively short evolution periods, complex ambient 
circumstances and gregarious nature. Two important approaches are systematic surveys 
and multi-wavelength studies towards individual sources to increase our knowledge about 
high-mass molecular cores that may harbor forming massive stars or mark the sites of future 
massive star formation. 

Previous surveys of high-mass star-forming regions focused on the sources associate d 
with ultra-compact (UC) Hll regions and their precursors (PUCHs) (IChurchwelll 120 02). 



Thes e works identified a number of high-mass protoste l lar objects (HM POs) ( IMolinari et al. 



1996. 2002: Sridharan et al. 2002: Bout her e1 al. 2002: \Yu el al. 20o3). However, the iden- 



tified objects usually have high luminosities (L IR > 1O 3 L ), indicating that most of them 
do not represent the earliest stage of massive star formation. By co mparing millimeter and 
mid-infrared (MIR) images of fields containing candidate HMPOs, ISridharan et al.l (120051 ) 
further identified a sample of potential high- mass starless cores (HMSCs), which may be 
the sites of future massive star formation. However, their MIR identification based on 
8.3 /im images from the Midcourse Space Experiment (MSX) was not sufficient to validate 
a genuine HMSC, because a heating accreting protostar may remain undetected up to 8 /im 
( IBeuther fc Steinackerl 120071 ) . The recently released high-resolution sensitive MIR and far- 
infrared (FIR) images obtained by the Galactic survey of the Spitzer Space Telescope could 
be used to verify these HMSC candidates. 

At the same time, some works towards individual sources indicated that massive cores 
at early stages might exist in the vicinity of evolved star- formi ng sites like UC Hll or Hll 
regions (IForbrich et al.l 12004 ; I Gar ay et al.l 12004 ; IWu et al.l 120051 ) . These works suggest that 
previously identified evolved sources may harbor objects at various evolution ary phases, in 



eludi ng HMSCs, high-mass cores harboring accreting protostars, and HMPOs (IBeuther et al. 



20071 ). One scenario is that the early-stage objects are stimulated by the star- forming ac- 
tivities in evolved regions. Another hypothesis is that they may form with their evolved 
companions during the fragmentation of parent clouds, but are restrained to give birth to 
stars in time by some physical supporting mechanisms. The third possible explanation is 
that the detected objects are just diffuse quiescent gas/dust clumps and will not form stars 
eventually. Probing the physical properties and circumstances of these objects may help to 
address the questions above. 



S87, cataloged as an optical H H nebula by ISharplessI (119591) . is a com plex star-forming 
region at a distance of ~ 2.3 kpc (Racine 19681 ; Crampton et al. 1978 ). It is associated 
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with a bright FIR s ource IRAS 19442+2427 and has been stu died by a numb er of authors. 
Henkel et al.l (119861 ) detected two 22 GHz water masers in S87. iBarsonyl (119891 ) studied it at 
radio, infrared, and optical wavelengths, suggesting the existence of a biconical outflow. A 
compact H H region was detected in centimeter radi o continuum, with an extended emission 



component (IBallv fc Predmore 



19831 ; lBarsonylll989l ). Two near-infrared (NIR) clusters were 



identified by I Chen et al.l (120031 ). labeled as S87E and S87W. The submillim eter (sub-mm) 
continuum emission of S87 exhibited an asy mmetric spatial configuration (IJenness et al. 



1995; Hunter et al. 2000; Mueller et a 



l ine map of CS J = 5-4 (IShirlev et al 



2002), which was also confirmed by the molecular 
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3). The previous works in ammonia (NH3) lines 



( IZinchenko et al.lll997l ; IStutzki et al.lll984l ) exposed two kinematically se parate comp o nents , 
which spatially overlap in the direction of S87E. The recent work of ISaito et al.l (120071 ) 
identified several gas clumps in the C 18 J = 1-0 map and proposed a hypothesis that 
S87E was formed by a cloud-cloud collision. All of the works mentioned above suggest that 
complex spatial and kinematic structures exist in S87, which may harbor objects at different 
evolutionary phases. The abundant data currently available from various wavelengths give 
us a great opportunity to perform a further comprehensive investigation towards S87. It 
may construct a consistent physical picture for this massive star-forming region and test the 
previously proposed hypotheses. 

In this article, we present a multi-wavelength study of S87, mainly based on the on- 
l ine archival data, our observations in m olecular lines, and two published observations 
( IZinchenko et al.lll997l ; iMueller et al.ll2002l ). We describe the used dataset and the observa- 
tional results of S87 in §|5]and[3J We concentrate on the spectral energy distribution (SED) 
analysis of the identified sub-mm clumps in §|H In §[5j we mainly discuss the stellar content, 
the star-forming activity and evolutionary stage of each sub-mm clump. The conclusions are 
summarized in §[6j 



2. DATA AND OBSERVATIONS 

2.1. Continuum Data 

All of the sub-mm/FIR/MIR continuum maps or images of S87 were obtained from data 
archives. 

The 850 and 450 fim sub-mm continuum data were retrieved from the James Clerk 
Maxwell Telescopd^ (JCMT) Science Archive, measured with the Submillimetre Common- 



1 The James Clerk Maxwell Telescope is operated by the Joint Astronomy Centre on behalf of the Sci- 
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User Bolometer Array (SCUBA) flHolland et all Il999f ) installed at JCMT. Two 850 and 
450 /um maps are available since S87 was observed twice in 2003. One observation was car- 
ried out in jiggle map mode on 2003 May 24th (JCMT program ID: M03AN23); the other 
was performed in Emerson II scan map mode on August 24th (M03BU45). The beamwidths 
of JCMT were 7". 5 (450 jttm) and 14"( 850/im). All of the retrieved data have been fully 
calibrated with the ORAC-DR pipeline ( Jenness et al. 20021 ) for flat-fielding, extinction cor- 
rection, sky noise removal, despiking, and removal of bad pixels, in the units of mJy beam -1 . 

The Spitzer MIR/FIR data were retrieved from the Spitzer Science CenteiB, includ- 
i ng the 3.6, 4.5, 5 .8, and 8.0 /im images measured with the Infrared Array Camera (IRAC) 
( IFazio et al.l 120041 ) and the 24 a nd 70 ittm images measured with the Multiband Imaging 
Photometer for Spitzer (MIPS) (TRieke et all l2004h . The IRAC and MIPS data are, re- 
spectively, from the G alactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) 
(IBeniamin et al.l 120031) an d the recently released MIPS Inner Galactic Plane Survey (MIPS- 
GAL) (jCarey et al.ll2005l ). All of them were calibrated by the Spitzer Science Center data 
processing pipel ines. In addition, we also retrieved the MIR images and point source catalog 
(PSC) of MSX ( lEgan et al.ll2003l ) from the Infrared Processing and Analysis Center (IPACjj 
for our study. 



2.2. Spectral Observation 

To investigate the molecular gas of S87, we mapped a region of 4' x 4' centered on 
IRAS 19442+2427 in the J = 1 -0 transitions of CO, 13 CO, C 18 and HCO+, with the 13.7m 
millimeter telescope of the Purple Mountain Observatory (PMO) in 2005 January and 2006 
May. A cooled SIS receiver was employed, and the system temperature T sys at the zenith 
was ~ 250 K (SSB). The backend included three acousto-optical spectrometers, which was 
able to measure the J = 1-0 transitions of CO, 13 CO, and C 18 simultaneously. All the 
observations were performed in position switch mode. The center reference coordinates are: 
R.A. (J2000) = 19 h 46 m 19 s .9, Dec. (J2000) = +24°35'24". The grid spacings of the CO and 
HCO + mapping observations were 60"and 30" respectively. The background positions were 
checked by single point observations before mapping. The pointing and tracking accuracy was 
better than 10". The obtained spectra were calibrated in the scale of antenna temperature T£ 



ence and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 

^http : / / ssc . spitzer . caltech . edu 



http : //www. ipac . caltech. edu 
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during the obse rvation, corrected for atmospheric and ohmic loss by the standard chopper 
wheel method (IKutner &: Ulichl Il98ll ). Table [1] summarizes the basic information about 
our observations, including: the transitions, the center rest frequencies u TeBt , the half-power 
beam widths (HPBWs), the bandwidths, the equivalent velocity resolutions (AV ms ), and the 
typical rms levels of measured spectra. All of the spectral data were transformed from the T£ 
to T m b scale with the main beam efficiencies before analysis. The uncertainty of brightness 
was estimated as 10%. The GILDAljjl so ftware package (CLASS/GREG) was used for the 
data reduction (IGuilloteau fc Lucasll2000l ). 



2.3. Other Archival Data 



We acquired the 350 /im continuum and NH 3 (J,K) = (1, 1) line maps of S87 through 
private communications with K. Young and I. Zinchenko. The 350 /im map was measured 
with the Sub-mm High Angular Resolution Camera (SHARC) installed at the Caltech Sub- 



mm Observatory (CSO) ( Mueller et al. 2002h . The NH3 (J, K) = (1, 1) line map was obtained 



with the Effelsberg 100 m telescope (IZinchenko et al.lll997l ). The technical details are sum- 
marized in the corresponding reference articles. 



3. RESULTS 

3.1. Sub-mm Maps 

Fig. [1] displays the 850 /im scan map (contours) and the Spitzer 8.0 /zm image (inverse 
grey-scale), in which the NIR clusters S87E and S87W are revealed as two bright MIR 
nebulae. The strongest peak of 850 /im is associated with S87E, and two other peaks exist 
to the northeast of it. We propose that these three 850 /im peaks are associated with three 
individual sub-mm clumps. They lie along an axis from southwest to northeast, hereafter 
labeled as SMM1, SMM2, and SMM3. 

We processed the sub-mm maps using the Richardson-Lucy (RL) iteration deconvolution 
algorithm to moderately enhance the spatial resolutions. As many other deconvolution 
solutions, this algorithm did not produce uncertainty information of the results. Therefore, 
we have to note that this process is not targeted to get the most "accurate" deconv olved 



maps. Our steps of deconvolutions are similar to those described by lSmith et al.l (120001 ) . The 
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850 and 450 /im beam patterns of SCUBA were constructed from the Uranus maps measured 
in 2003 August. A procedure in Starlink/KAPPA@ was used to perform the image processing 
tasks. We avoided the pixels at the edge of sub-mm maps during the iterations due to their 
low signal-to-noise level. 

The deconvolutions of the two 850 /im maps converged within 100 iterations and pro- 
duced acceptable enhanced maps without apparently artificial structures. However, for the 
450 /im maps, the procedure failed to converge within 150 iterations. Fig. [2] displays the de- 
convolved 850 /im maps and undeconvolved 450 /im maps. All of them have been converted 
into the units of mJy arcsec -2 . SMM3 is not covered in the jiggle maps (see Fig. [2fc and 
Fig. [2H) due to the limitation of the observational fields of view. SMM 1 is clearly elongated 
in the deconvolved 850 /im maps, and there are extended lobes to the west and south of its 
peak. SMM 2 is slightly elongated in the north-south direction. All of the three sub-mm 
clumps are revealed in a common envelop, suggesting that they may be associated although 
not necessarily in the same sky plane. 

To evaluate the CO J = 3 - 2 contribution to the 850 /im data, we examined our pre- 
vious observation of S87 in CO J = 3 - 2 at the Koln Observatory for Sub-mm Astronomy 
(KOSMA) 3m telescope. This observation was carried out for a CO multi-line survey of 
(UC) Hll regions and has not been published yet (Xue Sz Wu 2008, in preparation). After 
converting the integrated intensity of CO J = 3 - 2 to a flux density at 850 /im, we found that 
the contribution of CO J = 3 - 2 is less than the noise level of the 850 /im maps. We also 
evaluated the contribution of CO J = 6 - 5 to the 450 /im data by estimating its integrated 
intensity from CO J = 3 - 2 , under the assumption of local thermodynamic equilibrium 
(LTE) with an excitation temperature of 80 K. We found that its contribution is also small. 
Therefore, the effect of line contaminations can be ignored at 450 and 850 /im. 



3.2. Mid/Far-Infrared Images 

A luminous MIR point source is revealed at the position of the compact Hll region 
(see the IRAC images of Fig. [1] and Fig. [3]), which we henceforth label as MIRS A 
weaker MIR point source is found in the 3.6, 4.5, and 5.8 /tm images, ~ 8" to the south- 
west of MIRS 1. It is also detected by the 2MASS NIR All-Sky Survey and cataloged as 
2MASS 19461947+2435247. However, we did not find the NIR counterpart of MIRS 1 in 



E http : //www . j ach . hawaii . edu/ software/ starlink/ 

6 The recently released GLIMPSE Spring '05 Catalog names it SSTGLMC G60.8838-0.1282 
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2MASS images, indicating that MIRS 1 is highly obscured by the surrounding gas/dust en- 
velope at NIR wavelengths. In the zoomed 8 /jm image of Fig. [H two other point sources are 
found to the north of MIRS 1, which we label as MIRS 2 and MIRS 3. Strong diffuse MIR 
emission exists to the southeast of MIRS 1, coincident with the extended centimeter emission 



detected by iBarsonyl (119891 ). Faint diffuse MIR emission is detected in the northeast of the 
IRAC images, coincident with SMM3. SMM2 has no MIR counterpart in all of the IRAC 
bands. 

S87E and S87W saturate the 24 and 70 /im MIPS images. Five other sources are detected 
in the 24 /im band, which are coincident with the diffuse MIR emission in the IRAC bands. 
We label them as MIRS 4 to MIRS 8 (see Fig. EJ). Although the 24 and 70 /im images are 
saturated towards SMM1, it is still clear that the peaks of sub- mm and 24 /im emission are 
separate, which is also confirmed by a comparison with the MSX E band (21.3 /xm) image. 
SMM3 is associated with MIRS 4. However, its sub-mm peak and MIRS 4 are also slightly 
separate. There are only weak emission patches in the IRAC 5.8 and 8.0 /im bands towards 
SMM3, indicating that MIRS 4 is still embedded in its gas/dust cocoon. No 24 or 70 /im 
emission is detected towards SMM2, suggesting that it may be less evolved. 



3.3. Molecular Lines 



The ,7 = 1 — transitions of CO, 13 CO, and HCO + exhibit asymmetric line profiles 
(see the left panel of Fig. and two components are detected in C 18 ,7 = 1-0 . Since 
C 18 J — 1-0 is usually optically thin, we can rule out the possibility that the asymmet- 
ric li ne profile in the other transitions is caused by the self- abs orption in an i nf all enve- 



lope (IMvers et al. 



Zinchenko et al. 



1996; IWu et al.ll2005l ). Previous observations of IStutzki et al.l (119841 ) and 
(119971 ) also detected two separate compo nents in NHs (J, K) = ( 1,1) and 
NH 3 (J, K) = (2, 2) lines. The NH 3 (J, K) = (1, 1) spectra of Izinchenko et al.l Jl997h and our 
C 18 J = 1 - spectra at several positions are plotted in the right panel of Fig. HI These spec- 
tra further confirm that the broad lines of CO J = 1-0 , 13 CO J = 1 - , and HCO+ ,7=1-0 
consist of two components. 

We fitted our spectra at the reference position with Gaussian profiles. The results 
are displayed as the thin lines in Fig. HJ and the corresponding derived parameters are 
summarized in Table [21 including: the line center velocities, the fitted line widths, and 
the brightness temperatures. We estimated the beam-averaged column densities of C 18 
at the reference position using the standard LTE method. The excitation temperature of 
each component is assumed to be 35 K, in agreement with the estimation from CO J = 1 -0 
(assuming it is optically thick). The derived C 18 column densities are ~ 7.8 xlO 15 and 
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4.0xl0 15 cm 2 for the components at low and high velocities. 

Fig. [5^ is the HCO + J = 1 - position- velocity diagram along the northeast-southwest 
direction, which also exhibits two components. One is located at the reference position 
and associated with SMM1; the other extends from the reference position to the northeast, 
coincident with SMM2 and SMM3. We propose that these components arise from two 
clouds. Hereafter, they are named Cloud I and II, corresponding with the components at 
low and high velocity respectively. 

The integrated intensity maps of HCO + J = 1-0 and NH3 (J, K) = (1, 1) are also ex- 
hibited in Fig. O Two different integrated intervals are adopted, chosen to separate the 
emission from Cloud I and II. All of the presented intensity maps suggest: SMM2, SMM3, 
and the northeast part of SMM 1 may be associated with Cloud II; the main part of SMM 1 
is contributed by Cloud I. 



4. SED ANALYSIS 

4.1. Observational SEDs 

We extracted the 850 and 450 /im flux densities of each clump using a photometric 
procedure in the Starlink/GAIA software package. The measured results, as well as the 
positions and sizes of the adopted photometric apertures, are summarized in Table [3] We 
note that the uncertainties in Table [3] are just statistical errors (rms deviations derived 
from clean regions), and the estimation of the overall photometric uncertainties is difficult 
due to the limited information from the online data archives. However, a comparison among 
different observational modes and the previous similar observation may provide an evaluation 
of the accuracy of our results. 



Jenness et al.l (119951 ) observed S87 using the receiver UKT14 at JCMT in 1994. They 
detected two sources, which were coincident with SMM 1 and SMM 2 respectively. Our 
photometric results at 850 /im are in good agreement with theirs (see the last two columns 
of Table [3]), but the 450 /im results from SCUBA are systematically larger. Since UKT14 
is a single-element bolometer and its measurements may be affected by the change of sky 
conditions and other factors, we believe that the calibration of SCUBA data is more reliable. 
The photometric differences of the jiggle and scan maps are acceptable, less than 20% at 
450 /im. 



We examined the CSO map and the MIPS images to measure the flux densities of each 
clump at 24, 70 and 350 /im. Since SMM 1 saturates the 24 and 70 /im images, only lower 
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limits can be derived at these wavelengths. In addition, we checked the MSX PSC and found 
that the photometric apertures of SMM 1 and SMM 3 are coincident with the MSX point 
sources MSX6C G060.8828-00.1295 and MSX6C G060.9049-00.1275 respectively. Their flux 
densities are also adopted to construct the SEDs of SMM 1 and SMM 3. 

The measured flux densities are summarized in Table HJ extending from sub-mm to 
MIR. The average results of the scan and jiggle maps are adopted for 850 and 450 fim. 
Their differen ces are consid e red a s the uncertainties. The 350 fim uncertainties follow the 
description of iMueller et al.l (120021 ) and the uncertainties at 24 and 70 /im are the statistical 
errors. 



4.2. Isothermal Dust Model 



A simple isothermal gray-body d ust model is used to fit the observational SEDs. The 
details follow the method described by lSchnee et al.l (120071 ). In the adopted model, the mean 
weight of interstellar materials per hydrogen molecule \i is ~ 2.33. The dust opacity (mass 
absorption coefficient) K\ is dependent on the wavelength and can be described using the 
equation: 

1300 ftm.^ 



K\ — «130o(" 



A 



1 



where (3 is the dust opacity index and K1300 is the dust opacity at 1300 fim. Assuming a 
gas-to-dust ratio of 100, we adopt Ki.wi = 0.009 cm 2 g" 1 , which is derived from a gas/dust 



model with thin ice mantles (lOssenkopf &: Henningill994j ). 



We used a non-linear least-squares method (the Levenberg-Marquardt algorithm coded 
within IDL) to obtain the best-fit models for the observational SEDs. The physical properties 
of each sub-mm clump were obtained, including: the average dust temperature the dust 
opacity index /3, and the aperture-average column density of hydrogen molecules Nn 2 . In this 
fitting test, only the data upwards 70 //m were used. We assumed the 70 //m flux density of 
SMM 1 to be 4000 Jy, which was estimated from the interpolation of the IRAS flux densities 
subtracted with the potential contributions from SMM 2 and SMM 3. 

Fig. [6] displays the best-fit model SEDs for three clumps. We further calculated their 
clump masses and bolometric luminosities Lsed (by integrating the model SEDs over the 
range 1 fim< A < 2.0 mm). All of these derived re sults are summarize d in Table Addition- 
ally, we adopted a Monte Carlo method used by ISchnee et al.l (120071 ) to estimate the errors 
of the derived parameters that arise from the observational uncertainties. The 3a intervals 
are denoted as the superscripts and subscripts in Table [5j 
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The results in Table |5] show that the sub-mm clumps are all massive (110 — 220 M Q ). 
SMM 1 has a higher dust temperature, and its bolometric luminosity dominates in the whole 
region, implying the existence of strong internal heating source(s). The fitted dust opacity 
indices of three sub-mm clumps are slightly diffe rent (~ 1.3 — 1.8) and consistent with the 
typical values between 1 and 2 (IHill et al.1 120061 ). It must be noted that the derived Nn 2 
is directly affected by the adopted value of Kiwo- If we reduce K1300 by a factor of 2, Nn 2 
and the derived clump mass M, will increase by a factor of 2. However, the other derived 
parameters will not be affected by this change. 



4.3. Two-temperature Dust Model 

In Fig. [6j the best-fit models of SMM 1 and SMM 3 failed to describe the observational 
results below 70 /im. However, the model SED of SMM 2 can explain the absence of its 
MIR emission. To better characterize the excess MIR emission of SMM 1 and SMM 3, we 
performed another SED fitting test using a model with two dust components at different 
temperatures. In this fitting test, we adopted the observational data upwards 14.7 /im (ex- 
cluding 24 /im for SMM1). To reduce the fitting parameters, we assumed j3 is 1.5 and 1.3 
for each dust component of SMM 1 and SMM 3 respectively. The best-fit model SEDs are 
exhibited in Fig. [71 and the derived parameters of the warm and cool dust components are 
listed in Table El 

The two-temperature model fits the observational data very well above 12 /tin, which 
is consistent with the physical fact that there are warm dust around the internal heating 
sources and relatively cool dust envelopes surrounding the star-forming sites in SMM 1 and 
SMM 3. Although the IRAS 100 /tin flux density exceeds the model SED of SMM 1 (see 
Fig. [7J, we believe that the deviation is due to the large beam of IRAS . The results in 
Fig. [7] and Table [6] show that the warm components contribute little to the total masses and 
the flux densities at sub-mm wavelengths, but are required to explain the excess at MIR 
wavelengths. 

We tried to modify (3 to fit the emission below 12 /tm. However, no satisfying results 
were found. The emission in the MSX A and C bands does not follow the predication of 
gray-body models, suggesting these models are invalid at these wavelengths. Generally, two 
significant spectral features may exist at this MIR wavelength range. One is the emission 
of polycyclic aromatic hydrocarbons (PAHs), which is often detected towards compact Hll 
regions and photodissociation regions (PDRs). Previous studies have showed that the MSX 



A and C bands oft en contain PAH emission lines (IGhosh &: Ojhall2002l ; iKraemer et al.ll2003a 



Povich et al.1 120071 ). The other is the silicate feature, which has been predicted in the dust 
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model of lOssenkopf fc Henning (1994) an d demonstr ated to be important in the recent so- 
phisticated SED models (IRobitaille et al.ll2006l . 120071 ). This feature may be expected as the 
absorption at 9 .7 iim towards som e UC Hll regions, pro duced by the dust cocoons around 
center objects (IFaison et al.lll998l ). iPeeters et al.l (120021 ) have identified both of PAH and 
silicate features towards S87E in the previous ISO spectroscopy observation. The silicate 
absorption feature may be caused by Cloud II, which partly overlaps above the compact H II 
region. Since our gray-body models are focused to evaluate the overall properties of dust 
clumps, which are mainly constrained by the thermal emission from longer wavelengths, a 
detailed SED model explaining the PAH and silicate features is beyond our purpose. 



DISCUSSION 



5.1. Cloud-Cloud Collision 



All of the FIR/sub-mm images and molecular line maps exhibit the complex spatial 
and kinematic al structures of S 87. The recently published high-resolution C 18 J=l- 
observation (ISaito et al.l 120071 ) revealed several gas clumps at different velocities, which 
further confirms our i dentification of Clo ud I and II. However, are Cloud I and II really 
related to each other? ISaito et al.l (120071 ) proposed that the gas clumps at higher velocity 
might be on the near side along the line-of-sight because the observation of lChen et al.l (120031 ) 
detected many reddened sources in NIR there. They further reasoned that the clumps at 
low and high velocities were approaching and the NIR cluster S87E was possibly formed 
by a cloud-cloud collision. In the following, we verify the cloud-cloud collision model by a 
multi-wavelength comparison. 

Firstly, the 8.0 fim emission shows a sharp edge to the northeast of MIRS 1 (see Fig. [1]). 
This feature is probably caused by the large extinction at 8.0 //m because the sub-mm emis- 
sion is still strong. The position of the extinction patch is consistent with that of Cloud II, 
confirming that Cloud II is on the near side along the line-of-sight. Therefore, Cloud I and 
II are approaching. 

Next, we can infer from the intensity maps of Fig. [5] that the peak of SMM 1 and S87E 
are in the overlapping region of Cloud I and II. The MIR point sources in SMM 1 suggest 
that there is not only a formed NIR cluster, but also ongoing star-forming activities. The 
strong and continuous star-forming process is likely to be interpreted by the stimulation of 
a cloud-cloud collision rather than the spontaneous evolution of molecular clouds alone. 

Furthermore, the spatial c onfiguration of the compact Hll region and the associated 
extended centimeter emission (iBarsonyl Il989l ; iKurtz et al.l Il994f ) also supports the cloud- 
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cloud collision model. The champagne flow model of iKim fc Kool (120011 ) combined with 
clumpy structures of molecular clouds can explain the extended centimeter component which 
stretches to the southeast of the compact Hn region. If Cloud I and II are in contact, their 
contact plane will be along the northwest-southeast direction (see Fig. [5b and Fig. [5b). 
Consequently, the compact H II region will be better confined in the direction perpendicular 
to the contact plane and the champagne flow should be easier to spurt out in the southeast 
direction. If Cloud I and II are not in contact, the champagne flow will be more likely 
to splash in the direction perpendicular to the border of the parent cloud of the compact 
H II region. The observational result is consistent with the predication of the first scenario, 
supporting that the two clouds are colliding. Assuming that the two clouds have typical 
sizes R ~ lpc and a velocity separation 5v ~2kms _1 , the collision duration is at least 
R/5v ~ 5 x 10 5 yrs, comparable with the time scale forming a compact Hn region. 

The cloud-cloud collision is considered as an efficient mechanism to trigger star forma- 



tion. It may compress molecular gas and lead to local gravitational collapse (lLorenlll976 



Habe k Ohta 1992; Marinho et al. 



20011 ). However, its possibility is small in the diffuse 



molecular clouds (jElmegreenlll998l ) . Additionally, high velocity off- axis collisions cou ld be de- 
structive rather than lead to gravitational instabilities (jHausmanlll98ll ; lGilderull984f ) . There- 
fore, the fraction of star formation triggered by cloud-cloud collisions may be small in our 



Galaxy. All of the current evidence demonstrate s that S87 is a new examp 



collis i ons, and similar samples are still limited f|Loren 



1976 



Dickel et al. 



e of c 



1978 



oud-cloud 



Koo et al. 



19941 ; IVallee!ll995l : iBucklev fc Ward-Thomp^nlll996l : ISato et alihoool : lLoonev et al.ll2006h . 



5.2. Molecular Line Emission 



5.2.1. HCO+ J =1-0 



Our observation shows that the line profile of HCO + J =1-0 is similar to that of 
CO ,7 = 1-0 (see the left panel of Fig. 0J). Additionally, we found that both of CO J = 1-0 
and HCO + ,7 = 1-0 spectra show slight features of high- velocity (HV) gas when compared 
with C 18 ,7=1-0 , suggesting that HCO + extends in diffuse gas rather than simply con- 
centrates in the dense parts of gas clumps. Previous observational and theoretica l works 
have pointed out the abu ndance enhancement of HCO + in diffuse or shocked gas ([Turner 



1995bl ; iGirart et al.l 119991 ). which can explain our finding. 



However, the formation mechanism of the HV gas in S87 is unclear. We propose three 
different explanations: (i) the HV gas may arise from stellar outflows; (ii) it may be con- 
tributed by the high-pressure shocked material that is squirted out when the clouds collide; 
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(iii) or, it is from the non-impacting portions of the colliding clouds since they do n ot slow 



down to a common speed during the cloud-cloud collision. Although iBarsonyl (119891 ) identi- 
fied HV blue and red wings in her CO J = 1 - observation and proposed that the HV gas 
resulted from a biconical outflow with a wide opening angle viewed at large inclination, our 
identification of two individual clouds apparently rejects this model. High-resolution and 
sensitive observations are required to clarify the origin of the HV gas. 

Because both stellar outflows and cloud-cloud collisions can produce HV gas and broad 
non-Gaussian line profiles, it is possible that some observational results previously interpreted 
as bipolar outflows are caused by cloud-cloud collisions. However, since the possibility of 
cloud-cloud collisions is not high, similar cases like S87 should be rare. 



5.2.2. NH 3 (J,K) = (1,1) 



A feature of the NH 3 (J, K) = (1, 1) intensity maps is that the NH 3 emission tends to 
"evade" the luminous MIR sources. The NH 3 (J, K) = (1,1) peak of Cloud I is separate 
from MIRS 1 and the sub-mm peak of SMM 1. The NH 3 (J, K) = (1, 1) emission is absent to 
the southeast of MIRS 1, where the diffuse MIR emission is strong. The NH 3 ( J, K) — (1, 1) 
peak of Cloud I I is coincident wit h SM M 2, which ha s no M IR counterpart. In contrast, the 
observations of ISaito et al.l (120071 ) and IShirley et al.l (120031 ) showed that the C 18 J = 1-0 
and CS J = 5-4 emission is strong in SMM 1 and SMM 3. Since both of SMM 1 and SMM 3 
are dense clumps identified from sub-mm conti nuum, their rela tively weak NH 3 emission may 
be explained by the underabundance of NH 3 . iTurnerl (I1995al ) suggested that NH 3 could be 
destroyed by C + that dominates in PDRs. Howeve r, molecules like C* 8 Q are formed via C + 
and not affected by the photo-destruction process (IJansen et al.lll995l ). The diffuse 5.8 and 
8.0 fim emission near SMM 1 and SMM 3 is usually contributed by PAHs and interpreted as 
a tracer of PDRs. The existence of MIR emission there, as well as the strong C 18 J = 1 -0 
emission and the weak NH 3 emission, are consistent with the prediction of the chemical 
process proposed in previous works. 



5.2.3. virial States 



The line widths of molecular spectra are usually used to probe the kinematics of gas 
clumps. Since Cloud I and II can be well resolved in NH 3 lines, we estimate the virial masses 
of these two clouds in this section. 

We derived the line widths and brightness temperatures of NH 3 (J, K) = (1,1) at the NH 3 
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peaks of Cloud I and II, using the hyperfine structure fitting procedure of GILDAS/CLASS. 
The results are exhibited as thin lines in Fig. HI The angular diameters 8 Q i, s of Cloud I and 
II are calculated using the equation: 



9, 



obs 



(2) 



in which, Q is the measured angular area of each cloud. After that, we corrected the beam 
effect and estimated the intrinsic sizes of Cloud I and II following the equation: 



ft _ Y) V^cL ^mb 



(3) 



where R is the radius of the gas cloud in pc, D is the distance of S87, and 9 m b is the beamwidth 
of the NH3 ( J, K) — (1, 1) observation. Assuming that Cloud I and II are homogeneous 
spherical gas clouds with a density distribution p oc r~ a (a = 1.5) and neglecting the 
contributions from magnetic field and su rface pressure, the virial masses can be derived 
using the equation (IMacLaren et al.lll988l ): 



M vir = 126 ( 



5 -2a 



a: 



)RAV l 



FWHM) 



(4) 



in which, M v - n is the virial mass in M & and AVfwhm is the full width at half-maximum 
intensity (FWHM) of NH 3 (J, K) — (1, 1) in km s -1 . 

All the measured and derived parameters of Cloud I and II are listed in Table in- 
cluding: the positions of NH 3 peaks, the angular and intrinsic sizes, the line widths at NH 3 
peaks, and the derived virial masses of two clouds. The total virial mass of Cloud I and 
II is ~430M©, which is much smaller than the previous estimati on (~1080M Q ) obtained 
with the total line width of two components (IZinchenko et al.lll997l ) but comparable with the 
mass estimated from the SED fitting (~460M Q , from the isothermal dust model). However, 
we note that the above comparison of the masses estimated from different approaches can 
be affected by the adopted dust opacity and the assumption of a. Although a variation of a 
is not likely t o cause much change in viria l masses, the dust opacity may change by at least 
a factor of 2 (jOssenkopf fc Henning!ll994l ). which leads to a large uncertainty in the masses 
estimated from SEDs. 



5.3. Stellar Contents of SMM 1 and SMM 3 

The position of MIRS 1 is consistent with that of the compact Hll region, within the 
astrometric error (1.5"), indicating that MIRS 1 is the exciting massive (proto)star. We 
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examined the high-resolution centimeter map of iBarsonyl (119891 ) and found that neither 
MIRS 2 nor MIRS 3 shows compact radio continuum emission. The possible explanation is 
that MIRS 2 and MIRS 3 are less evolved compared with MIRS 1 or t hey are not massive 
enough to ionize their surroundings and to excite compact H II regions. iHenkel et al.l (119861 ) 
detected a strong water maser in SMM1, which is often considered to be associated with 
HMPOs. The velocity range of this water maser is 21 — 25kms _1 , in good agreement with 
the systematic velocity of the molecular clouds. All the evidence mentioned above supports 
that SMM 1 is a high-mass star-forming site that harbors massive forming stars or cluster. 

The Lym an continuum radi ation from massive stars mainly escapes in the form of free- 
free emission. iKurtz et al.l (11994 ) estimated that the Lyman continuum photon flux required 
to keep the entire region of S87E ionized was 3.2 x 10 46 photons s" 1 , which correspo nds to that 
of a B O. 5 ZAMS star. The bolometric luminosity of such stars is ~ 3.0 x 10 4 L Q (ICrowther 
20051 ). slightly smaller than that of SMM 1 (~ 3.7 x 1O 4 L , from the two-component model). 
The extra luminosity of SMM 1 may come from the relatively weak MIR sources near MIRS 1, 
which can not be traced by the free-free emission. 

SMM 3 contains the bright 24 fim source MIRS 4. The ratio of the luminosities from its 
cool and warm components is ~ 2.1, lower that of SMM1. Its bolometric luminosity is ~ 
740 Lq, also lower than that of SMM1, which indicates that SMM 3 is more likely to be an 
intermediate-mass star-forming site. 



5.4. Physical Properties of SMM 2 



No MIR point source or diffuse emission below 70 jum is detected towards SMM 2. Since 
only a cold dust component can describe its observational SED, strong internal heating 
sources are not likely to exist in SMM 2. 



arises 



Henkel et al.l (119861 ) detected a weak 22 GHz water m aser near SMM 2, which usually 
from the dense circumstellar disks arou nd protostars (Park &: Choil2007 ) or originates 



in outflows from the birth of a massive star ( van Dishoeck fc Blake! 1998 ). Since this water 
maser is in the velocity range 8 — 15kms _1 , significantly different from that of the molecular 
clouds, we favor the second explanation for its origin. We notice that this water maser is on a 
sub- mm emission ridge connecting the peaks of SMM 1 and SMM 2 rather than near the peak 
of SMM 2, and its position uncertainty is large when compared with sub-mm observations. 
Therefore, we doubt that this weak maser is produced by the intrinsic factors of SMM 2. For 
instance, the potential outflows from massive protostars of SMM 1 may shock the ambient 
molecular gas of SMM2 and produce a weak water maser at the rear side of SMM2. This 
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scenario is consistent with the lower velocity of the water maser. Therefore, we believe that 
the existence of this water maser does not necessarily contradict SMM 2's physical properties 
derived from the SED and MIR image analyses. Based on the information available, we 
support that SMM 2 is probably a HMSC that may form massive stars or intermediate star 
clusters eventually. 

6. CONCLUSIONS 

We have carried out a multi-wavelength study of the massive star-forming region S87. 
The main results are summarized as follows. 

1. We identified three sub-mm clumps in S87, labeled as SMM1, SMM 2, and SMM 3. 
They are estimated to have masses of 210, 140, and 110 M , with average dust temperatures 
of 41, 21, and 24 K respectively (from the isothermal gray-body model). 

2. We examined molecular line maps from our observations and compared them with 
previous results of other authors. We concluded that the star-forming activities in SMM 1 
are stimulated by a cloud-cloud collision. 

3. We found that HCO + can trace diffuse gas and NH 3 may be destructed by chemical 
processes in the region harboring MIR sources or exhibiting strong diffuse MIR emission. 

4. We calculated the virial masses of the two colliding clouds, which are in good agree- 
ment with those estimated from SEDs. 

5. The stellar contents and star- forming activities of sub-mm clumps are identified. 
Their SEDs reveal that these clumps are at various evolutionary stages. SMM 1 and SMM 3 
are high-mass and intermediate-mass star-forming regions respectively. SMM2 is massive 
and cold, has no MIR counterpart, which is probably a HMSC. All of these results expose 
that the star formation in S87 is at multiple phases. 

We are grateful to K. Young and I. Zinchenko for their sharing of the 350 /im and 
NH 3 ( J, if) = (1, 1) maps. We would like to thank the staff at the Qinghai Station of 
PMO for their assistance during the observations and Weilai Gu for her help to obtain 
the HCO + J = 1-0 data. We acknowledge the anonymous referee for his/her careful reading 
and helpful suggestions. Qifeng Yin and Sophia Day are also thanked for their help on the 
manuscript. This research was funded by the Grant 10733030 and 10128306 of NSFC. It 
employed the facilities of the Canadian Astronomy Data Center operated by the National 
Research Council of Canada with the support of the Canadian Space Agency, and it is also 
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Table 1. Observation Parameters 



Transition 




^rest 


HPBW 


Bandwidth 


AV xes 


lcr rms a 






(GHz) 


(") 


(MHz) 


(km s _1 ) 


(K chan. -1 ) 


CO J = 1 - 





115.271204 


46 


145 


0.37 


0.18 


13 CO J = l 





110.201353 


47 


43 


0.11 


0.14 


C 18 J = l- 





109.782182 


48 


43 


0.12 


0.14 


HCO+ J=l- 





89.188521 


58 


43 


0.26 


0.11 



a typical value in the scale of T£ 



Table 2. Parameters of Molecular Lines 



Transition 3 






AV 










(km s- 1 ) 


(km s _1 ) 


(K) 


13 CO J = 


1- 


o 1 


21.51 (0.08) 


1.93 (0.04) 


16.71 (0.38) 


C 18 J = 


1- 


o 1 


21.48 (0.08) 


1.69 (0.19) 


2.39 (0.39) 


13 CO J = 


1- 


o 2 


23.61 (0.04) 


2.60 (0.06) 


15.04 (0.38) 


C 18 J = 


1- 


o 2 


23.78 (0.16) 


1.83 (0.42) 


1.14 (0.39) 


HCO+ J = 


1- 


o b 


22.30 (0.04) 3.10 (0.16) 


4.54 (0.11) 



*The error levels are from Gaussian fitting. 

a The superscripts indicate the different components. 

Since the two components can not be well resolved, the 
derived parameters are from single Gaussian fitting. 



Table 3. Photometric Results of the Sub-mm Clumps 



Object 




S 


#ap 


850 ^m a 


450 ^m a 


850 nm h 


450 fj,m h 


850 fim c 


450 /im° 




J2000 


J2000 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


SMM1 


19 46 19.8 


f 24 35 32 


60 


17.7±0.1 


133±10 


16.8±0.1 


167±1 


17 


110 


SMM2 


19 46 22.3 


f 24 36 01 


30 


5.8±0.1 


32±6 


5.3±0.1 


45±1 


6.5 


18 


SMM3 


19 46 23.2 


f 24 36 28 


30 


4.1±0.1 


25±6 











a flux densities extracted from the scan map (M03BU45) 
b flux densities extracted from the jiggle map (M03AN23) 
c flux densities from lJenness et al.l (|1995f l 



-22 - 



Table 4. Flux Densities of the Sub-mm Clumps 



Object 


850 Aim 


450 /im 


350 /an 


70 /im 


24/nm 


21.3 nm a 


14.7 /zm a 


12.1 ^m a 


8.3 /zm a 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


SMM1 


17.3±0.9 


150±34 


256±46 


> 1500 b 


> 24 b 


225±13 


43.3±2.3 


31.7±1.6 


19.6±0.8 


SMM2 


5.5±0.6 


39±13 


79±14 


49±3 


< 3 C 










SMM3 


4.1±0.6 


25±6 


40±8 


36±3 


9.1±3.0 


5.1±0.4 


0.6±0.1 


l.lzbO.l 


1.3±0.1 



a These values are derived from MSX PSC. We adopted the flux densities of 
MSX6C G060.8828-00.1295 and MSX6C G060.9049-00.1275 for SMM1 and SMM3 respec- 
tively. SMM2 is not detected in all of the four MSX bands. 

b Only the lower limits are available due to the saturation of the Spitzer 24 and 70 /xm 
images. We assume the 70/itm flux density of SMM 1 to be 4000±200 Jy when performing the 
SED fitting. 

c no detection 



Table 5. Best-fit Parameters for the SEDs (I) 



Object 


T d (3 
(K) 


Nn 2 a 
(10 22 cm- 2 


M 
) (Af ) 


^SED b 

(L Q ) 


SMM1 
SMM 2 
SMM 3 


4itg i.5i2:2 
2il? i-8±8:| 

24t 2 2 1.3±8-J 


q 9+0.9 
°' z -0.5 

8.418:2 
6.4±8:g 


210122 

140+18 

lioisg 


31000 
1100 
600 



a aperture-average column density. 

b Thc luminosities are calculated from the integra- 
tion under the model SED curves over the range 
1 /im< A < 2.0 mm. 
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Table 6. Best-fit Parameters for the SEDs (II) 



Object 


Id 


P 


iv H2 a 


M 


T b 
^SED 


— component 


(K) 




(10 20 cm" 2 ) 


(M Q ) 


(L e ) 


SMM 1— cool 


40 


1.5 


320 


210 


30000 


— warm 


92 


1.5 


0.7 


0.5 


7200 


SMM 3— cool 


23 


1.3 


690 


112 


500 


— warm 


82 


1.3 


0.4 


0.1 


240 



a aperture-average column density. 

b The luminosities are calculated from the integra- 
tion under the model SED curves over the range 
1 /zm< A < 2.0 mm. 



Table 7. Parameters of virial Mass Estimatation 



Object 


Vlsr 




AS a 


Ay 


#obs 


R 


vir 




(km s" 1 ) 


(") 


(") 


(km s- 1 ) 


(") 


(pc) 


(Mq) 


Cloud I 


21.01(0.06) 








1.22(0.15) 


126 


0.67 


170 


Cloud II 


23.84(0.02) 


40 


40 


1.60(0.05) 


116 


0.60 


260 



*The error levels are from the spectral fitting. 

a The positions of intensity peaks are listed as offsets 
from the reference coordinates: R.A. (J2000) = 19 h 46 m 19 s .9, 
Dec. (J2000) = +24°35'24". 



The virial mass estimated by IZinchenko et al.l (|1997D is 1080 L 



which was derived from the total width of two components in the 
NH 3 (J, K) = (1, 1) spectra. 



-^4:36:00.0 



30.0 
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35:00.0 




O Henkdetal.(1986) 
□ Barsony (1989) 



24.0 



22.0 



_l I 1_ 



19:46:20.0 18.0 
, I i i i 1_ 



S8 



« 



O JCMT/S CUBA 850 




Fig. 1. — JCMT/SCUBA 850 /xm continuum emission (contours) overlaid on the 
Spitzer /IRAC 8.0yum image. The contour levels increase from 0.6 to 5.0 Jy beam -1 , in 
a step of 0.4 Jy beam -1 (4a). The inverse grey-scale 8.0 /im image is in the unit of MJy 
sr -1 . The c oordinate sys t em is J2000. The two small open circles mark the water masers 



detected by lHenkel et al.l (119861 ). The square denotes the compact Hn region, and the open 
ellipse indicate s the position and rough size of the associated extended centimeter emission 
(lBarsonylll989l ). The open circle at the bottom shows the beamwidth of JCMT at 850 /im. 
The sub-mm clumps SMM1, SMM2, and SMM3, as well as the NIR clusters S87E and 
S87W, are labeled. Three 8.0 /im point sources, named as MIRS 1, MIRS2, and MIRS3 in 
§ 13.21 are marked by the arrows. 
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Fig. 2. — Contour plots of the sub-mm continuum maps, overlaid on their inverse grey-scale 
images. The beamwidth is indicated in the bottom-right corner of each panel, (a) The 
deconvolved 850 /im scan map (M03BU45), after 98 iterations. The contours increase from 
1.8 to 29 mJ arcsec -2 in a step of 3.1 mJ arcsec -2 . The three open circles over the emission 
peaks denote the photometric apertures (see §HJ). (b) The undeconvolved 450 /im scan map 
(M03BU45). The contours are: 40, 120, 220, and 380 mJ arcsec" 2 . (c) The deconvolved 
850 yum jiggle map (M03AN23), after 51 iterations. The contours increase from 2 to 26 mJ 
arcsec -2 in a step of 3 mJ arcsec -2 . (d) The undeconvolved 450 /im jiggle map (M03AN23). 
The contours are: 60, 120, 200, 280, 360, 440, and 500 mJ arcsec -2 . 
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Fig. 3. — The Spitzer MIR/FIR images at six wavelengths. The sizes and positions of six 
panels are the same. The absolute coordinates (J2000) are exhibited in the 8.0 /im panel. 
The white areas in the emission regions of 24 and 70 /im images are due to the saturation of 
the detectors. The MIR emission is also plotted as contours in the 3.6, 4.5, 5.8, and 8.0 fim 
panels. Their contour levels are: 1%, 2%, 3%, 5%, 10%, 16%, 28%, 50%, and 90% of the 
peak intensities. The crosses in each panel denote the 850 /im peaks. The square indicates 
the compact H II region. The contour plot of the 24 /im panel is the 850 /im continuum 
emission, adopted from Fig. [2h. MIRS 4 to MIRS 8 are also marked in the 24 fim panel. The 
small open circles in the bottom-right corner of 24 and 70 //m panels denote the resolutions 
of these observations. 
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Fig. A.— Left panel: The spectra of CO J =1-0 , 13 CO J =1-0 , C 18 7=1- 
, and HCO+ J =1-0 at the reference coordinates: R.A. (J2000) = 19 h 46 m 19 s .9, 
Dec. (J2000) = +24°35'24". The vertical dashed lines define the border of the C 18 7 = 1-0 
emission above 3a level. Right panel: The C 18 Q 7 = 1-0 sp ectra from our observation and 
the NH 3 (7, if) = (1,1) spectra from IZinchenko et al.l (11997T ) . The spatial positions of the 
spectra are denoted in the top-left corners of small plots. The thin lines on the spectra 
are the fitted results. We performed double Gaussian fitting for the 13 CO 7 = 1-0 and 
C 18 7 = 1-0 spectra. For the NH 3 (7, if) = (1,1) transition, we used the NH 3 hyperfine 
structure fitting procedure of GILDAS/CLASS. 
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Fig. 5. — In each plot, the dashed contour denotes the 3a level. The center reference coordi- 
nates are: R.A. (J2000) = 19 h 46 m 19 s .9, Dec. ( J2000) = +24°35'24". (a): The HCO+ J = 1-0 
position-velocity diagram along the northeast-southwest direction. The contour levels in- 
crease from 3a to 21 a by 3a. The vertical dashed line indicates the border of two integrated 
intervals (see in § !3.3p . (b) and (c): The integrated intensity maps of HCO + J = 1-0 . 
The integrated interval is showed in the top- right corner of each plot. The contour levels 
increase from 40% to 100% by 15% (5<r) of the peak value. The square denotes the com- 
pact H II region, and the open ellipse indicate s the position and rough size of the associated 
extended centimeter emission f lBarsonyl Il989l ). The crosses denote the 850 /im peaks, and 
the beamwidth is indicated in the bottom-right corner of each plot. The dashed line is the 
location of the potential contact plane mentioned in § 15.11 The position of S87W is marked 
by triangles, (d) and (e): The integrated intensity maps of NH 3 (J, K) = (1, 1) . The c ontou r 
levels increase from 3a in a step of 3a. The spectra data are from IZinchenko et al.l (119971 ). 
The positions of MIRS 1 and MIRS4 are marked as triangles. 
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Fig. 6.— From (a) to (c), the best-fit SEDs (solid lines) of SMM1, SMM 2, and SMM3 
are exhibited. The filled triangles denote the flux densities of IRAS 19442+2427. The open 
triangles denote the flux densities extracted from MSX PSC. The open squares indicate the 
flux densities derived from the Spitzer /MIPS images. The small filled squares on the model 
SEDs are the data points used for the SED fitting. The errorbar of each data point is plotted 
if available. 
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Fig. 7.— The panel (a) and (b) exhibit the best-fit SEDs of SMM 1 and SMM 3 (solid lines). 
The two components of model SEDs are also displayed (dot lines) respectively. The filled 
triangles denote the flux densities of IRAS 19442+2427. The open triangles denote the flux 
densities extracted from MSX PSC. The open squares indicate the flux densities derived 
from the Spitzer /MIPS images. The small filled squares on the model SEDs are the data 
points used for the SED fitting. The errorbar of each data point is plotted if available. 



